More common errors in finding exact solutions 
of nonlinear differential equations. I 

Roman 0. POPOVYCHt* and Olena 0. VANEEVA* 

t Fakultat fur Mathematik, Universitdt Wien, Nordbergstrafle 15, A-1090 Wien, Austria 
$ Institute of Mathematics of NAS of Ukraine, 3 Tereshchenkivska Str., Kyiv-4, Ukraine 
E-mail: rop@imath.kiev.ua, vaneeva@imath.kiev.ua 

In the recent paper by Kudryashov [Commun. Nonlinear Sci. Numer. Simulat., 2009, 14, 
3507-3529] seven common errors in finding exact solutions of nonlinear differential equa- 
Q\ ' tions were listed and discussed in detail. We indicate two more common errors concerning 

the similarity (equivalence with respect to point transformations) and linearizability of dif- 
ferential equations and then discuss the first of them. Classes of generalized KdV and 
mKdV equations with variable coefficients are used in order to clarify our conclusions. We 
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investigate admissible point transformations in classes of generalized KdV equations, ob- 

Qtain the necessary and sufficient conditions of similarity of such equations to the standard 
KdV and mKdV equations and carried out the exhaustive group classification of a class of 
variable-coefficient KdV equations. Then a number of recent papers on such equations are 
commented using the above results. It is shown that exact solutions were constructed in 
these papers only for equations which are reduced by point transformations to the standard 
KdV and mKdV equations. Therefore, exact solutions of such equations can be obtained 
from known solutions of the standard KdV and mKdV equations in an easier way than 
by direct solving. The same statement is true for other equations which are equivalent to 
well-known equations with respect to point transformations. 
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Finding exact (resp. closed-form, resp. explicit) solutions plays a significant role in investigation 
of physically important partial differential equations, especially, nonlinear ones. There exist a 
number of famous named solutions in different fields of physics. First such solutions were con- 
structed and intensively studied for models describing flows of fluids (e.g., the Navier-Stokes and 
Euler equations and the boundary layer equation) since it is difficult to neglect nonlinear effects 
arising in fluid dynamics. These were the Poiseuille flow, the Burgers vortex, Jeffery-Hamel solu- 
tion for the flow in a wedge shaped region, the von Karman solution for flow over a single rotating 
disc etc. (see, e.g., the extensive reviews by Berker and by Pukhnachev [21 [25]). Simultaneously 
with penetrating the nonlinear paradigm into other fields of science, the construction of exact 
solutions became more and more topical and involved nonlinear diffusion and wave equations 
and systems of such equations, nonlinear Schrodinger equations, general relativity equations etc. 
Some equations have only a few known exact solutions (e.g, the Ablowitz-Zeppetella traveling 
wave solutions of the Fisher equation), for other equations wide families of exact solutions pa- 
rameterized with a number of arbitrary constants and functions have been obtained (integrable 
and linearizable equations and equations possessing large Lie symmetry groups). Last decades 
the scientific activity in this field permanently grows. 

Recently jTT] Kudryashov has made an excellent review and listed seven common errors 
in finding exact solutions of nonlinear differential equations, which appeared in the modern 
mathematical and physical literature. The above errors relate to correctness, generality and 
representation of such solutions and methods of their constructions. In spite of its carefulness, 
the list by Kudryashov does not exhaust all existing common errors in this subject. We indicate 
two more common errors, continuing numbering of [IT] for convenience of further references. 
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These errors usually arise on the very initial step of consideration, namely, under the choice of 
differential equations for finding exact solutions and related investigations. 

Eighth error. Exact solutions are often constructed with no relation to equivalence of differ- 
ential equations with respect to point (resp. contact, resp. potential etc.) transformations. 

Ninth error. Exact solutions of linearizable differential equations are sometimes found with 
no relation to the linearizability property. 

Mathematical investigations are always carried out up to equivalence relations defined in 
sets of similar objects. Equivalent objects differ only in properties which are inessential for 
certain consideration and, therefore, can be identified in certain sense. This rule is extended to 
mathematical models, arising in applications. 

A usual equivalence used for differential equations is that generated by locally nondegener- 
ate point transformations and called similarity [19]. Similar systems of differential equations 
have similar local properties and similar related local objects (local solutions, symmetries, local 
conservation laws etc.). 

The purpose of this letter is to discuss the eighth error. In fact this error can be considered 
as a whole family of different inaccuracies concerning similarity of differential equations. We 
distinguish, e.g., the following kinds of inaccuracies. 

1. Any differential equation from the class under consideration is similar to the same clas- 
sical differential equation for which wide multiparametric families of exact solutions were 
constructed earlier and possibilities for finding new solutions inequivalent to known ones 
look as problematic. 

2. Only differential equations from a subclass of the class under consideration are similar 
to classical ones and exact solutions are obtained only for such equations. Sometimes 
constraints singled out the subclass of equations with found particular solutions from the 
whole class are explicitly presented, sometimes they are implicitly implied or even not 
indicated in any way. 

3. The class of differential equations under consideration can be mapped by a family of 
point transformations parameterized with arbitrary elements of the class to another class 
of simpler structure (resp. with less number of arbitrary elements, resp. with simpler 
transformational properties etc.) but this possibility is missed. 

4. Similarity of differential equations is used for finding exact solutions but the usage is 
incomplete (resp. improper, resp. incorrect). 

There are rather a lot of published papers in which exact solutions were constructed for 
equations similar to well-known and/or simpler equations without using the similarity. Choos- 
ing examples for commenting, we restricted ourselves with quite recent papers published in 
certain journals and devoted to considerations of equations related to the KdV and mKdV 
equations. Our choice is justified by the fact that both the equations are the most known and 
well- investigated nonlinear equations of mathematical physics. They are connected via the Miura 
transformation. Since these equations are integrable and have nice symmetry properties, wide 
families of their exact solutions had been found using different powerful methods (the inverse 
scattering method, Darboux and Backhand transformations, the Hirota bilinear method etc.) 
and collected in a number of papers, textbooks and handbooks (see, e.g., [U HU H7J ESI EZ] and 
references therein). At the same time, in the commented papers only very particular results on 
similar equations were obtained. 
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The illustration of the above list of inaccuracies with generalized KdV and mKdV equations 
is additionally justified by their own importance in different fields of physics (see, e.g., references 
in commented papers). Moreover, the existence of similarity transformations between equations 
which are related to the KdV and mKdV equations is a well-known fact from the late 1960's. 
A simple point transformation connecting the Gardner and mKdV equations was already pre- 
sented in [T7]. A more nontrivial point transformation between the standard and cylindrical 
KdV equations was found in [13] (see also [7]). Then the similarity of certain generalizations 
of KdV and mKdV equations with coefficients depending on t or (t, x) to these equations was 
investigated in detail (see, e.g., [H EJ [8] and references therein). Later the similarity ar- 
guments were permanently used under reviewing and refereeing papers similar to commented 
ones, cf. [15] . Lie and generalized symmetries of such equations were also investigated (see, e.g., 

!3 123130!). 

Our paper is organized as follows: In the next section we construct a hierarchy of normalized 
classes of third order (1 + l)-dimensional evolution equations, which is related to the examples 
commented. This gives the complete description of admissible point transformations within such 
classes. For two wide subclasses of the variable-coefficient KdV and mKdV equations, jointly 
covering equations from almost all the examples, we derive the necessary and sufficient condition 
of similarity of such equations to the standard KdV and mKdV equations, respectively. Exhaus- 
tive group analysis of a normalized class of variable-coefficient KdV equations is carried out in 
Section O Section [J] is partitioned into two subsections. In subsection 14.11 we collect examples 
when the whole considered classes of equations are similar to the KdV or mKdV equations, 
exemplifying the first of the inaccuracies listed above. Subsection 14.21 includes examples on 
particular similarity, which illustrate the second and third inaccuracies. Different errors under 
the construction of similarity transformations among generalized KdV and mKdV equations are 
classified and illustrated by examples in Section In the conclusion other possible applications 
of similarity transformations are discussed. 



2 Admissible point transformations in classes of generalized 
KdV equations 

Following [20j , we start from the general class of third order evolution equations and construct a 
hierarchy of nested normalized subclasses of this class, which consist of different generalizations 
of the KdV and mKdV equations. In this way we describe the entire sets of admissible point 
transformations of these subclasses. At first we briefly define necessary notions, considering 
only the case of single differential equations in one unknown function for simplicity. Rigorous 
definitions for the general case of systems of differential equations are given, e.g., in [20J. 

Let Cg be a differential equation L(y,u^,9(y,u^)) = for the unknown function u of 
n independent variables y = (yi, . . . ,y n )- Here U(p\ denotes the set of all the derivatives of u 
with respect to y of order not greater than p, including u as the derivative of order zero, and 
L is a fixed function depending on y, UU)\ and 9. The tuple 6 of k arbitrary (parametric) 
functions 9 1 (y,u^), . . . , 9 k (y,u^) runs through the set S of solutions of the auxiliary system 
S(y> u (p),9(q)(yi u (p))) = of differential equations with respect to 9. In this system y and 
play the role of independent variables and 9r q \ stands for the set of all the partial derivatives 
of 9 of order not greater than q with respect to the variables y and u^y Usually the set S is 
additionally constrained by the condition Yi(y,Uf p \,9r q \(y,urp\)) ^ with another differential 
function E. In what follows we call the functions 9 arbitrary elements. Also, by C\s we denote 
the class of equations Cg with the arbitrary elements 9 running through S. The set T(£|s) 
of triples (9,0,ip), where 9,9 £ S and ip is a point transformation from Cg to Cg if such a 
transformation exists, is called the set of admissible transformations in the class C\s- (This is 
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a formalization of the notion of allowed [30] or form-preserving [TO] transformations.) Roughly 
speaking, the (usual) equivalence group G~ of the class C\s consists of the point transformations 
in the space of variables and arbitrary elements, which are projectable on the variable space and 
preserve the whole class C\s- The class C\s is normalized if any admissible transformation in 
this class is induced by a transformation from G~. 

Consider the general class £ 3 of third order evolution equations. They have the form 

Uf — H(t, X, U, U x , U X xi ^xxx)i 

where H Uxxx ^ 0. To find the set of admissible transformations of the class £ 3 and its complete 
equivalence group including both discrete and continuous equivalence transformations, we apply 
the direct method. The calculations are simplified with taking into account the well-known 
fact that the expression for t in any point (and, even, contact) transformation connecting two 
(1 + l)-dimensional evolution equations depends only on t |10|. I14j . Thus, a point transformation 
maps an equation from <? 3 to an equation from the same class if and only if it has the form 

t = T(t), x = X(t, x,u), u = U (t, x, u). (1) 

The functions T, X and U have to satisfy the nondegeneracy assumption TfA ^ 0, where 
A = X X U U — X U U X . The equivalence group Gq of £ 3 consists of the transformations 

, A U t D x X - X t D x U 

t = T(t), x = X(t,x,u), u = U(t,x,u), H = jrj^H + f , (2) 

where T, X and U run through the corresponding sets of smooth functions satisfying the above 
nondegeneracy assumption and D x denotes the operator of total differentiation with respect to 
the variable x, D x = d x + u x d u + ut x d Ut + u xx d Ux + • • • . Therefore, the class £ 3 is normalized. 

The subclass £q.i singled out from £ 3 by the constraint H UxxxUxxx = has the same equivalence 
group Gq and, therefore, is normalized. The same claim is true for the subclass £q 2 °f ^oi' 
associated with the additional constraint H UxxxUxx = 0. 

Setting one more constraint H UxxxUx = leads to the normalized subclass £ 3 X of £q 2 whose 
equivalence group G± is a proper subgroup of Gq consisting of the transformations ([2]) with 
X u = 0. The normalized subclass £f 2 singled out from £f A by the constraint H UxxxU = has 
the same equivalence group G^. 

The equivalence group G 2 of the normalized subclass £ f nested in £ 3 2 and associated with 
the additional constraint H UxxxX = is properly contained in G~. Its elements additionally 
satisfy the condition X xx = 0. 

Another possibility is to impose the constraint H Uxx = within the subclass £qi- The 
corresponding subclass £f is normalized and possesses the equivalence group G^ formed by the 
transformations ([2]) for which X u = 0, U uu = and 2U UX X X = U U X XX . 

Then the subclass £\ = £ 2 n £f is normalized and its equivalence group is G% = G 2 n . 
Integrating the total set of the constraints imposed on the arbitrary element H and the conditions 
obtained for elements of G% gives that the equations from the class £ f have the form 

u t + g{t)u xxx = F(t,x,u,u x ), 

where g ^ 0. The equivalence group G± consists of the transformations (we present only 
components corresponding to the equation variables) 

t = a(t), x = (3(t)x + j(t), u = 9(t)u+$(t,x), (3) 

where a, 0, 7, 9 and <3? are arbitrary smooth functions of their arguments, a{59 ^ 0. 
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Further we consider two more special subclasses of E\ which are directly related to examples 
commented in Sections H] and 

The first subclass is formed by the variable-coefficient KdV equations 

u t + f{t)uu x + g{t)u xxx + h(t)u + {p{t) + q{t)x)u x + k(t)x + l(t) = 0, (4) 

where all the parameters are arbitrary smooth functions of t, fg ^ 0. This subclass is normalized. 
Its equivalence group (in terms of the arbitrary element H) is singled out from by the 
condition & xx = 0. Hence the components of transformations corresponding to the equation 
variables have the simple general form 

t = a(t), x = f3(t)x + -y{t), u = 6(t)u + <p(t)x + ip(t), (5) 

where a, (3, 7, 6, ip and ip run through the set of smooth functions of t, a/39 7^ 0. The arbitrary 
elements of @ are transformed as follows 

/ = — 3/, 9 = —9, h = — \h--f 

ato at at \ u 

~ 1 ( Vf±&\ ~ 1 (r , 7^~/^ , , Pt 
q = — 9-77/ + ^" ' P = — \PP ~ 79 H n / + 7t - 7 



a< V $ P / &t \ 6 ' P 

i = ^(«-** + «) + £/-<p. + <4). 

f=^(«-^(»*-* + «) + T / -'" + ''?)-*(''"? / "^)- w -* 

Any equation from class (jll) can be reduced by point transformations to the form ([8]) with g = 1 
and forms (|10p and (|15p . The necessary and sufficient condition of similarity of such equations 
to the standard KdV equation is 



St = 2gs — ?>qs H k, where s := 1 — „ . (6) 

9 9 19 

The second subclass consists of the variable-coefficient mKdV equations 

u t + f{t)u 2 u x + g(t)u xxx + h(t)u + (p(t) + gCtJxJiia + fc(t)iw a + Z(t) = 0, (7) 

where all the parameters are arbitrary smooth functions of t, fg 7^ 0. This subclass is also 
normalized. Its equivalence group (in terms of the arbitrary element H) is singled out from G% 
by the condition § x = 0, i.e., the components of transformations corresponding to the equation 
variables are of the form ([5]) with cp = 0, where a, P, 7, 8 and i/j run through the set of smooth 
functions of t, aP9 7^ 0. The arbitrary elements of ([7]) are transformed by the formulas 

f P f ~ P 3 I 1 ( h 0t 
/=— a>f> 9 = —9, h= — [h- — 
att) at at \ a 

+ -r 1 P = — \Pp ~ 79 + - P^k + 7 t - 7- 



a t V P J ott V # P 

i= £(*- 2 ?')- r-£(«-*-*+*2 

Five of the arbitrary elements can be gauged to simple constant values. For example, it is 
possible to set g = 1 and h = p = q = l = 0. An equation of form (|7|) is similar to the standard 
mKdV equation if and only if 

2/i-2g = 4"-, 2lf = k t + kh-k^. 

f 9 f 
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In Section 0] we present different examples on similarity of equations from classes @ and (JT]) 
to the standard KdV and mKdV equations, respectively. 

Note that hierarchies of normalized classes in the same superclass, which are different from 
the constructed hierarchy, can be easily obtained according as purposes of investigation. 



3 Group analysis of a class of variable-coefficient KdV equations 

We study a class of variable-coefficient KdV equations in detail since this creates a basis for the 
consideration of four examples in Section 14.21 The general form of these equations is 

ut + f(t)uu x + g{t)u xxx = 0, (8) 

where / and g are arbitrary (smooth) functions of t, fg ^ 0. 

As shown in Section [21 for the general values of the parameter-functions (arbitrary elements) 
/ and g equation ([8]) is not equivalent to the standard KdV equation up to point transformations 
but at least one of the parameters (/ or g) can be set equal to 1 using a point transformation. 
(The equations of the form (JHJ with g = 1 are called the transitional KdV equations [3].) Thus, 
after the transformation 

t = J f(t)dt, x = x, u = u (9) 

equation ([8]) takes the same form with f(t) = 1 and g{t) = g(t)/f(t). Therefore, without lost of 
generality we can consider the class of equations 

u t + uu x + g(t)u xxx = 0, (10) 

where g is an arbitrary (smooth) nonvanishing function of t, since this class is the image of 
class (JHJ) under the mapping generated by the family of point transformations ([9]). (See [201 129] 
for related definitions.) 

We carried out the exhaustive group classification of class (|10p . Note that Lie symmetries 
and admissible ("allowed") transformations of the wider class of equations having the form 

u t + f(t,x)uu x + g(t,x)u xxx = (11) 

were investigated in [30]. The separate consideration of subclass (|10p is additionally justified 
by the fact that it has nicer transformational properties than the superclasses ([8]) and, espe- 
cially, (nj>. 

The equivalence group G~ of class (llOj) consists of the transformations 
~ at + b _ e2X + e\t + eo 



ct + d ' ct + d 



(12) 



eAct + d)u — C2cx — cqc + e\d e2 g 

e ct + de 

where a, b, c, d, eo, &\ and are arbitrary constants, e = ad — be ^ 0, and without loss of 
generality we can assume that e = ±1. The group G~ is isomorphic to the matrix group 

&2 7^ 0, ad — be = ±1 

Equations from class (|10p are similar only if they are G~ -equivalent. Moreover, all admissible 
transformations in this class are generated by transformations from G~, i.e., the class (jlOp is 
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normalized in the usual sense. (The initial class dH) is normalized only with respect to the 
extended generalized equivalence group and class (fTTj) possesses no normalization properties.) 
This implies the following claim: An equation of form (llOj) is similar to the KdV equation if 
and only if gtt = 0. Any transformation realizing the similarity belongs to G~. Therefore, an 
equation of form ([8]) is reduced to the KdV equation by a point transformation if and only if 

g(t) = f(t) L J f(t) dt + co) , (13) 

where Co and c\ are constants, (cq,ci) / (0,0) [4J, that well agrees with Theorem 3 of [30] (cf. 
also equation ©). Equation (|13p coincides with the constraint on arbitrary elements of the 
equations from class ([8|) which have the Painleve property [9]. 

Denote by A 9 the maximal Lie invariance algebra of the equation of form (|10p with the fixed 
value g of the arbitrary element. It follows from the infinitesimal invariance criterion |X8|. [T9] 
that the coefficients of any vector field Q = r(t, x, u)c\ + £ (i, x, u)d x + r](t, x, u)d u from A 9 satisfy 
the system of determining equations 

r x = t u = Cu = r]uu = Cxx = i] xu = 0, j] = (£ x - T t )u + r] t + r] x u = 0, (14) 
Tgt = (3&C - r t )g. 

Integrating equations (|14p gives the following expressions for the coefficients of Q: 

T = C 2 t 2 + C\t + C , 

£ = (c 2 t + c 3 )x + c 4 t + c 5 , 

rj = (-c 2 t + c 3 - c\)u + c 2 x + c 4 , 

involving only arbitrary constants Co, . . . , C5, i.e. for any g 7^ the algebra A 9 is finite dimen- 
sional. Then the last determining equation leads to the unique classifying condition 

(c 2 t 2 + cit + co)g t = (c 2 t + 3c 3 - ci)g. 

The group classification of class (fTUj) is equivalent to the integration of the classifying condition 
up to the G~-equivalence. The transformation (112j) acts on the coefficients of the equation 
(at 2 + j3t + j)g t = (at + S)g in the following way: 

a = a 2 a + ac(3 + c 2 7, 

= 2aba + (ad + bc)/3 + 2cdj, 

7 = b 2 a + bd(3 + d 2 j, 

5 = aba + bc(3 + cd7 + eS. 

The kernel ^4 n = Dg^oA 9 of the maximal Lie invariance algebras of equations from class (|10p 
is A n = (d x , td x + All G~-inequivalent cases of Lie symmetry extension are exhausted by 
the following: 

1. g = exp ( / .f^tf, dt] , (a, 5) + (0, 0), (a, (3, 7) # (0, 0, 0): 
\J at 2 + /ft + 7 / 

,49 = ^ + a M) 3( Q t2 + /3t + 7 )d t + (3 Q f + S + (3)xd x + ((-3at + 5- 2[5)u + 3ax)a u ). 

Up to G~ -equivalence this case is partitioned into the three inequivalent subcases: 

g = e 1 : A 9 = (d x , td x + d u , 3d t + xd x + ud u ) 
(a = j3 = 0, S/j is scaled to 1); 



7 



g = tP,H> 1/2, (J, + 1: A3 = (d x , td x + 8 U , 3td t + (/i + + (/z - 2)ud n ) 

(a = 7 = 0, M := 8/0); 

g = e SaKbaat y/W+l: A° = (d x , td x + d u , 3(t 2 + l)d t + (3t + <5)x^ + ((-3* + + 3x)d M ) 
(a = 7 = 1, = 0); 

2. 5 = 1: A 9 = (d x , td x + d u , d t , 3td t + x9 x - 2ud u ). 

The presented group classification gives all inequivalent values of g for which the classical 
method of Lie reduction can be effectively used. 
The class flSJ) can be also mapped to the class 

u~ t + uu x + u xxx + h(t)u = (15) 

by the family of point transformations 

t = J g{t) dt, x = x, u = —u. 

The arbitrary element /i of the mapped class is expressed via the arbitrary elements / and g in 
the following way: 

f _ f(t)g t (t)-ft(t)g(t) 

[ ) ~ f(t)(g(t)) 2 • 

Thus, the equation (|10p with g = t is then mapped to the cylindrical KdV equation (h = (2i) _1 ) 
whose similarity to the standard KdV equation is known for a long time [13]. Analogously, the 
value g = e* corresponds to the spherical KdV equation (/t = which is not integrable. 



4 Examples on similarity of equations and their solutions 

As shown in Section [21 to establish similarity of equations in all the presented examples it 
is sufficient to apply point transformations of the simple general form ([5]). In some of the 
examples presented we change the notations of arbitrary elements in equations in order to unify 
the consideration. Thus, the arbitrary elements /, g, h, k, I, p, q and s run through the set of 
smooth functions of t satisfying certain conditions of nonvanishing, which are explicitly indicated 
in the examples. We also explicitly indicate necessary conditions of nonsingularity for arbitrary 
elements even if they were not given in the corresponding commented papers. We should like to 
emphasize that correctness of solutions given in commented papers is not discussed here. Other 
common errors of different kinds described in [11] were made in the most of commented papers. 

4.1 Similarity of entire classes 

In each example collected in this section, all equations of the class under consideration are 
similar to a single well-studied equation (either the KdV or mKdV ones). This means that 
the arbitrary elements (constant or functional parameters) introduced in order to generalize the 
classical equations to such classes are completely needless. 

Example 1. In [C28] the so-called "generalized (G'/G)-expansion method" was applied to 
finding exact solutions of "the mKdV equations with variable coefficients" having the form 

u t = K (t)(u xxx - 6u 2 u x ) + 4JCi(t)u x - h(t)(u + xu x ), (16) 

where Kq, K\ and h are arbitrary (smooth) functions of t, Kq ^ 0. The same class of equations 
was investigated, e.g., in [C2j . |C5| and |C16j using so-called "the variable-coefficient generalized 



8 



projected Ricatti equation expansion method", "generalized Jacobi elliptic functions expansion 
method" and "sub-ODE method" , respectively. 

At the same time, any equation of form (|16p is equivalent to the standard mKdV equation 
with respect to a point transformation. Indeed, consider the point transformation 

u 

t = a(t), x = P(t)x + 1 (t), u = e—, (17) 

where a, (3 and 7 are smooth functions of t, which should be additionally specified, at/3 7^ and 
e £ { — 1; 1}. It maps equation (fTl)j) to the equation 

a t ui + -j{u + (3xu x ) + "ftUx = K /3 3 (uxxx - 6u 2 tii) + 4Kx(3u x - h(u + (3xu z ). 

If, e.g., e = 1 and the functions a, and 7 satisfy the system 

a t = -K /3 3 , (3 t = -h(3, it = 4ifi/3, 

the transformed equation is nothing but the standard mKdV equation u^—6u 2 u x + u X x$ = 0. In 
other words, the function u = u(t, x) satisfies equation (|16p if and only if there exists a solution 
u = u(t, x) of the standard mKdV equation such that 

u = (3u(a, j3x + 7), 

where a = - f K e~ 3 f hdt dt, /3 = e -f hdt , 7 = 4 / K x e~ J hdt dt. 

Note additionally that the usual equivalence group [20], H9] of class (fl~6l) consists of the trans- 
formations of form (|17p prolonged to the arbitrary elements Kq, K\ and h by the formulas 



at at V 44\ /3/y at \ [3 



Moreover, if two fixed equations from the class (|16p are connected via a point transformation 
then this transformation has the form (|17p . This means that the class (|16j) is a normalized 
class of differential equations which is a single orbit of its equivalence group, generated from the 
standard mKdV equation. 

Example 2. Some traveling wave solutions of the "compound/combined KdV-mKdV equation" 

u t + (a + (3u)uu x + 'jUxxx = 0, (18) 

where a, (3 and 7 are real constants, (3^ ^ 0, were constructed in \C3\ IC71 IC131 IC2H IC24} IC29|, 
IC32] . In fact, this equation is called the Gardner equation (a should be scaled to a standard 
value) and is obviously similar to the mKdV equation Uf+£u 2 u x + Uxxx = 0, where e = sign(/?7), 
with respect to the point transformation 



a 2 



t = jt, x = x + —t, 



a 



which is well known for a long time [T7]. Therefore, each solution of equation (|18p is represented 
in the form 



u{t, x) 



2 

a \ a 



where u is a solution of the mKdV equation, and for any solution u of the mKdV equation this 
representation gives a solution of (fT8|) . 
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The other close class of equations 



u t + u x + au u x + u xxx = 0, (19) 

where a runs through the set of real nonvanishing constants, was considered in |C25] . Only 
specific traveling wave solutions were found using the so-called "Exp-function method". The 
equations of form (|19p were groundlessly called the modified Benjamin-Bona-Mahony equations. 
In fact the Benjamin-Bona-Mahony equation ut + u x + uu x — u xx t = has no, at least, direct 
relation with class (I19p . At the same time, any equation of form (119j) is reduced by the trivial 
point transformation 



t = t, x = x — t, u = y\a\u 

to the mKdV equation u^ + eu 2 u x + u xxx = 0, where e = sign a. Equation (fT9j) with a = 1 was 
also investigated using "extended F-expansion method" in |C9j . 

Example 3. The authors of [CT1] apply "generalized expansion method" to find exact solutions 
of generalized KdV equations with variable coefficients, which have the form 

u t + g(t)(6uu x + u xxx ) + Qf{t)g(t)u = x(f t (t) + I2g{t)f 2 (t)) + M(t) (20) 



with g 7^ 0. The whole class (|20j) is mapped to the KdV equation + Guu x + u xxx = by the 
family of point transformations 

gj 3 dt, x. = 72 — 6 J g y 3 (3 dt, u = ^ 0, 

where 7 = e~ 6 -l fgdt and (3 = J M-f~ 2 dt. This means that the function u = u(t, x) satisfies an 
equation of form (|20p if and only if it is represented via a solution u of the KdV equation by the 
expression 



An auto-Backlund transformation and exact solutions of equations from class (|20p with M = 
were considered in [C18j . All related results also can be reproduced from corresponding results 
for the KdV equation using the above point transformations, where M = 0. 

The subclass of the equations (|20p with g = 1 and M = arose in |Clj . where symmetry 
properties of such equations were studied. It was also mentioned in |Clj that obtained results 
can be extended to the equations of the more general form 

u t + 6uu x + u xxx + 6f(t)u = x(f t (t) + 12/ 2 (t)) + ht(t) + I2f(t)h(t). (21) 



The family of point transformations mapping class (|21h to the KdV equation consists of the 
transformations 

/o [ u — fx — h 
7 dt, x = jx — 6 j hjdt, u = 2 ) 

where 7 = e -6 -^ $ dt . 

Example 4. Using the so-called "Jacobi elliptic function expansion method" , some new soliton- 
like solutions were obtained in [C31| for the KdV-like equations with variable coefficients, which 
have the general form 

u t + f{t){xu x + 2u) + g(t)(u xxx + auu x ) + h(t)u x = 0. (22) 
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Here a is an arbitrary constant and ag ^ 0. Any equation from class (|22p is similar to the 
standard KdV equation ui + uu x + u xxx = with respect to the point transformation 

t = a(t), x = f3(t)x + 7(t), u = (gj^2 u > ( 23 ) 



where the parameter- functions a, [3 and 7 are expressed via arbitrary elements of class (|22 [) in 
the following way: 



(3 = e -ff dt , a = Jgp 3 dt, j = -Jh(3dt. 



4.2 Particular similarity 

Although the classes of evolution equations from examples of this section are not similar in 
whole to classical equations, they also involve needless arbitrary elements. Canceling such 
arbitrary elements with point transformations leads to essential simplification of calculations. 
Moreover, in all the papers commented here results were in fact derived only for specific values 
of arbitrary elements for which the corresponding equations are similar to either the KdV or 
mKdV equations. 

Example 5. Class ([HJ) of KdV-like equations with coefficients depending on t was considered 
in \C26\ IC30] using the so-called "Exp-function method" . (See also references in [C26J IC30| on 
a number of papers where some particular results were obtained for the same class in similar 
frameworks.) Auto-Backlund transformations and some analytical solutions for equations from 
class (|H1) were considered in |C4l IC6} IC17j . In fact, certain results are obtained in the above 
papers only for the particular case g/f = const of (fT3|) which is obviously reduced to the 
KdV equation [4]. The corresponding point transformations combined with the known wide 
families of exact solutions of the KdV equation give at once much more knowledge about exact 
solutions of equations from class ([8]) than that presented in |C26} IC30| and previous papers on 
the subject. In [C20] Lax pairs were constructed only for equations whose arbitrary elements 
satisfy condition (fl~3j) . This is why the results of [C20| are obvious consequences of classical 
results for the KdV equation. See the discussion on Lax pairs in the conclusion. 

Example 6. In [ClOj "generalized variable coefficients KdV equations" of the form 

ut + f{t)uu x + g(t)u xxx + l(t) = (24) 

with fgj^O were considered using a "new generalized expansion method" . The coefficient I can 
be made equal to zero by a point transformation. More precisely, any equation from class (|24f) 
is mapped to the equation from class dH) with the same values of the arbitrary elements / and g 
by the transformation 

t = t, x = x + r y(t), u = u + ip(t), 

where 74 = ftp and tpt = Since class (f24"|) is a preimage of class (|8|) and, therefore, class (fl~0|) 
with respect to mapping generated by families of point transformations, it is needless to construct 
exact solutions of equations of form (|24p with nonzero values of I. 

Moreover, exact solutions of (124p were constructed in [CIO] only in the case 

3gft 2 ~ 3ff t g t - fgftt + f 2 9tt = 

when the corresponding equation (|24h is obviously similar to the KdV equation in view of that 
the above equation in / and g is equivalent to condition (fT3|) . 



11 



Example 7. At first sight the class of "KdV equations with variable coefficients" 

u t + f(t)uu x + g{t)u xxx + h(t)u + (p(i) + q(t)x)u x = (25) 

considered in |C91 IC19[ [C27J looks quite complicated since it involves the five (!) arbitrary 
smooth functions /, g, h, p and q of t, where fg ^ 0. 

The form (|25p was partially simplified in |C27] in the beginning of consideration. Namely, the 
coefficient h was gauged to zero by a point transformation. In fact only one arbitrary element 
among /, g, h, p and q is essential. The other ones can be set equal to standard fixed values. 
Indeed, class (|25|) is mapped to class ([8]) by the family of point transformations 



t = t, x = (3{t)x + 7(t), u = 9(t)u, 

where 6 = e^ hdt , (3 = e~f qdt and 7 = — f p(3dt. The corresponding value of the imaged 
arbitrary elements are / = and g = g/3 3 . The class ([8]) can be then transformed, e.g., to 

class (fTUj) or class (fT5]h 

Exact solutions were constructed in [C191 IC27| only for equations on form (125j) with arbitrary 
elements constrained by the condition 

fe^ h - 2q)dt /g = const 

which is written down in terms of the imaged arbitrary elements as f /g = const. In |C9| the 
constraint h = 2q was additionally imposed. Hence all such equations are similar to the KdV 
equation (cf. Example [5]) and the construction of their exact solutions is needless. 

Example 8. The class of systems of the form 

Ut + h(t)uu x + s(t)vv x + g(t)u xxx = 0, 

(2b) 

v t + f{t)uv x + g{t)v xxx = 

generalizing the Hirota-Satsuma system were studied in [C33] by using so-called "F-expansion 
method". Here f,g,h and s are arbitrary smooth functions of t satisfying the conditions 

fgs / 0, g/f = const, f -h = a 2 s, 

where a is a constant. The parameter-function g could be at once set equal to 1 by a simple 
transformation only of t (namely, t = J gdt), then the transformed element / would be an 
arbitrary nonvanishing constant. 

The solution presented in |C33| satisfies the additional condition v = —au under which the 
system (|26j) is reduced, in view of the condition / — h = a 2 s, to the single equation (jSJ) which 
is similar to the KdV equation since g/f = const (cf. Example [5]). As a result, the solutions 
of (|26p found in [C33] in fact can be easily derived from simplest traveling wave solutions of the 
KdV equation. 

Example 9. In [C23| some exact solutions were constructed for variable coefficient mKdV 
equations of the general form 

u t + f(t)u 2 u x + g{t)u xxx + p(t)u x = 0, (27) 

where fgy^O and additionally the condition gj f = const should be satisfied. At the same time, 
any equation of form (|27p . whose coefficients are constrained by this condition, is mapped by 
the point transformation 



t = / 9(t) dt, x = x- p(t) dt, u = a/I/7 



ig\u, 
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to the mKdV equation + eu 2 u x + u xxx = 0, where e = sign(fg). Without additional suppo- 
sition, an equation from class (|27p is reduced by point transformations to an equation from the 
same class in which, e.g., p = and / = 1. Another possibility is to gauge the coefficients p and 
g to the values and 1, respectively. 
The more general class of equations 

u t + f(t)u 2 u x + k(t)uu x + g(t)u xxx + p(t)u x + h(t)u = (28) 

with fg^O was considered in |C22| within the framework of Hirota's approach. Only the equa- 
tions of form (|28p whose arbitrary elements satisfy the additional constraint e~ 2 f hdt f j g = const 
were represented in bilinear forms. One- and multi-solitary- wave solutions of such equations 
were found after imposing the one more additional constraint e~f hdt k/g = const. (Under the 
construction of multi-solitary- wave solutions this constraint was used implicitly.) 

Exact solutions for the equations from class (|28|) were also constructed in |C19j using the 
so-called "extended tanh method". The solutions also were obtained only for those equations 
which coefficients satisfy the system of both above constraints. 

Any equation from class (|28|) can be mapped to an equation from the same class, in which 
p = 0, h = and g = 1. The corresponding transformation of the variables and other arbitrary 
elements is 

-2 fhdt - fhdt 

p(t)dt, u = ef hdt u, f = /, k = k. 

9 9 

Then the above constraints take the form / = const and k = const, respectively. If the arbitrary 
element / is constant, it can be set to ±1 by scaling u. If we additionally have k = const, 
this parameter can be made to vanish by a simple transformation (cf. Example [2]). Therefore, 
bilinear forms were given in [C22J only for equations similar to the equations 

u t + eu 2 u x + k(t)uu x + u xxx = 0. 

Moreover, in [C191 IC22| exact solutions were found only for equations similar to the mKdV 
equation! 

Note that the above constraints as well as canonical forms of the corresponding equations up 
to point transformations had already been presented, e.g., in |C8| (cf. also Section [2|). 

Example 10. Lie symmetries and similarity solutions of generalized KdV equations of the 
general form 

u t + u n u x + g(t)u xxx + h(t)u = 0, (29) 

where n £ Z + and g ^ 0, were investigated in [C12j . The parameter- function h can be set equal 
to zero by the point transformation 

' /*"'"• 

where 6 = ef hdt and the transformed value of the arbitrary element g is g(t) = g(t)(9(t)) n . 
This means that fixing the arbitrary element h cannot lead to cases of equations (|29p with 
special symmetry properties. Therefore, results from [C12] on Lie symmetries of equations from 
class ([29]) are at least incomplete. 



g(t) dt, 
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5 Inaccuracies concerning similarity transformations 



In spite of that the similarity transformations in classes of variable-coefficients equations gener- 
alizing classical integrable ones are well known [U El [201 EO] , in a recently published series of 
papers (see, e.g., |C81IC14llC15| and references therein) methods of the construction of similarity 
transformations between different generalized KdV (resp. mKdV, resp. nonlinear Schrodinger 
etc.) equations were presented as quite original. Although these results are not really new, the 
above presentations of them contain certain common inaccuracies which can be summed up in 
the following points: 

• Restrictions on involving arbitrary elements for which found similarity transformations are 
correct are not indicated in an explicit form. 

• Equations are reduced to a few canonical forms which are similar to each other. 

• Wide sets of transformations reducing equations to canonical forms are constructed in 
cases when it is enough to look for only the simplest one. At the same time, no proof on 
the completeness of the sets is presented. 

• Cases of reducibility of equations to canonical forms are partitioned into different subsets 
although such partition is needless. 

• Similarity transformations are looked out only in special sets of points transformations 
without rigorous arguments. 

For each of the above points there is a simple justification. Thus, equations whose arbitrary 
elements satisfy different conditions of, e.g., vanishing or nonvanishing can be reduced to non- 
similar canonical forms. If two canonical forms C\ and £2 are similar with respect to a point 
transformation Tq then an equation L is reduced to C\ by a point transformation T\ if and 
only if the same is true for £2 with the transformation TqTi. A point transformation T maps 
an equation C to an equation £ if and only if it can be represented in the form T = TqS = 
S'Tq, where Tq is a fixed point transformation from C to S and S' are point symmetry 
transformations of C and respectively. Given a class of differential equations and a canonical 
form which some of these equations are similar to, it is better to present a unified condition which 
is singled out the entire subclass of equations reducible to the canonical form. Most physically 
important classes of differential equations admit only fiber-preserving point transformations 
satisfying different simple conditions (see, e.g., Section [2] and [101 l2"0] ). 

Example 11. Class (|28p was considered in |C8| without indicating explicit conditions on arbi- 
trary elements. It was shown that an equation from this class possesses the Panleve property if 
and only if the corresponding values of the arbitrary elements satisfy the additional constraints 
e -2fhdtf/ g _ cons t anc [ e - I hdt qjg = const. Lax pairs, Backhand transformations and simplest 
solutions of such equations were found. Then point transformations which reduce the equations 
possessing the Panleve property to the Gardner and mKdV equations were constructed. 

Such consideration has a few inaccuracies. The general form (I28D should be simplified from 
the very beginning (cf. Example [9]). Depending on the sign of f/g, an equation of form (|28p with 
the Panleve property is similar either the KdV equation or the mKdV equation with the same 
sign of the nonlinearity. This is why at least the condition / 7^ should be explicitly imposed 
for the mKdV equation to be a canonical form. The equations with the Panleve property 
were reduced also to the Gardner equation which is obviously similar to the mKdV one (cf. 
Example [2]) and, therefore, contains, as a canonical form, a needless arbitrary element. The 
unnecessary parameters appeared also in constructed transformations. 

Analogous remarks can be made on |C14| . where equations of the form with k = q = 
and fg^O were considered. Under proper additional constraints on arbitrary elements, such 
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equations were simultaneously reduced to two canonical forms (the standard and cylindrical KdV 
equations) which are similar. Transformations were a priori assumed to be foliation-preserving 
although the necessity of this assumption can be easily proved. The similarity transformations 
constructed should be factorized with respect to point symmetry transformations. 

In [C19| it was incorrectly claimed that any equation from class (125p can be reduced to the 
cylindrical KdV equation. 

The similarity of equations of the form 

u t + F(t, x)uu x + G(t, x)u xxx + H(t, x)u + P(t, x)u x + L(t, x) = (30) 

with L = to the KdV equation was considered in |C15| . Only foliation-preserving point trans- 
formations which are linear in u were chosen for establishing the similarity, and the choice was 
not justified. This approach was groundlessly called "the direct method". The system of de- 
termining equations for components of the transformations was presented in an inaccurate way 
with no simplification and integrated only under a strong additional constraint. As a result, all 
equations the similarity of which to the KdV equation was really studied belong to the class . 
At the same time, it can be easily proved (cf. Section [2|) that the class (I30p is normalized. Its 
equivalence group consists of the transformations whose components corresponding to the equa- 
tion variables have the form ([3]). The necessary and sufficient conditions of the similarity to the 
KdV equation (which are analogous to ©) and the corresponding point transformations can be 
simply derived for both the entire class (I30j) and its subclass singled out by the constraint L = 0. 
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6 Conclusion 

We have discussed the common error in finding exact solutions of nonlinear differential equations, 
which concerns the similarity of equations. The consideration is exemplified with a number of 
papers on variable-coefficient KdV and mKdV equations. There also exist flows of similar papers 
on the variable-coefficient sin-Gordon and Schrodinger equations, different generalizations of 
diffusion equations etc., which should be additionally analyzed. Note that a part of such papers 
concerning nonlinear Schrodinger equations have been already examined in |12t I20|. 
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The similarity can be applied not only for generating new solutions from known ones and 
simplifying calculations. It also gives different criteria for testing final results. Thus, if an arbi- 
trary element of a class of differential equations can be neglected using a point transformation, 
it should play no crucial role in conditions singling out, from this class, subclasses with special 
symmetry or transformational properties or special cases in finding exact solutions. 

The approach based on similarity is easily extended to different local objects and properties 
related to differential equations, e.g., Lie and point symmetries |201 1191 129]- conservation laws 
and potential symmetries [211 E21 E3] , reduction operators (i.e., nonclassical symmetries) pHEH], 
Backlund transformations etc. For example, given a "spectral problem" V: ip x = Mip, ipt = Vip 
associated with an equation C (with no boundary conditions), a point transformation T: t = 
T(t), x = X(t, x), u = U(t, x, u) with T t X x lI u ^ maps V to the system V: ip x = Uip, tp^ = Vi[), 
where U = X x ~ l U and V = Tt~ 1 (V — X x XtU), whose compatibility condition C is similar 
to C with respect to T. At the same time, even similar equations can demand a separate 
investigation, e.g., of related boundary problems and asymptotic properties [16]. 

The next paper of the series will be devoted to the ninth common error. We plan to exem- 
plify the discussion on finding exact solutions of linearizable differential equations with papers 
on equations from the Burgers hierarchy, their different generalizations and multidimensional 
versions. 
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